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Abstract Using 64 ms count data for long gamma-ray bursts ( Tgo > 2.6 
s), we analyze the quantity, relative spectral lag (RSL), which is defined 
as T3i/FWHM(^i), where T31 is the spectral lag between energy bands 1 
and 3, and FWHAI/^i-f denotes full width at half maximum of the pulse in 
energy channel 1. To get insights into features of the RSLs, we investigate in 
detail all the correlations between them and other parameters with a sub- 
sample including nine long bursts. The general cross-correlation technique 
is adopted to measure the lags between two different energy bands. We can 
derive the below conclusions. Firstly, the distribution of RSLs is normal 
and concentrates on around the value of 0.1. Secondly, the RSLs are weakly 
correlated with FWHM, asymmetry, peak flux (Fp), peak energy {Ep) and 
spectral indexes {a and /?), while they are uncorrelated with rai, hardness- 
ratio {HR31) and peak time The final but important discovery is that 
redshift (z) and peak luminosity (Lp) are strongly correlated with the RSLs 
which can be measured easily and directly. We find that the RSL is a good 
redshift and peak luminosity estimator. 
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1 INTRODUCTION 



The temporal profiles of gamma-ray bursts (GRBs) generally exhibit very complex and 
variable characteristics due to overlapping between adjacent pulses (Norris et al. 1996; 
Quilligan et al. 2002). So far many investigations on the analysis of their light-curves 
especially the pulses have been made. For example, the properties of pulses such as 
widths, amplitudes, area of pulses and time intervals between them together with number 
of pulses per burst had been studied by several authors (e.g. McBreen et al. 1994, 2001, 
2003; Li et al. 1996; Hurley et al. 1998; Nakar & Piran 2002; Qin et al. 2005). 

In addition, some investigations associated with spectra have also been made (e.g. 
Kouveliotou et al. 1993; Hurley et al. 1992; Ghirlanda et al. 2004a). In particular, the 
spectral lag between variation signals in different energy bands not only reflects the 
features of spectrum evolution but also exhibits the properties of light-curve. Many re- 
searches regarding this variable had been done from many distinct aspects (see, e.g. 
Norris et al. 2000, 2001, 2005; Gupta et al. 2002; Kocevski & Liang 2003; Daignc & 
Mochkovitch 2003; Schaefer 2004; Li et al. 2004; Chen ct al. 2005). It is interestingly 
found that unlike short GRBs the spectral lags of most long GRBs arc larger than zero 
and concentrate on the short end of the lag distribution, near 100 ms ( Band 1997; Norris 
et al, 2001). 

Concerning the redshift (or luminosity) indictors with GRBs, previous investigations 
have offered us some significant paradigms in the case of light-curves, for instance, the 
relations of luminosity-lag (Norris et al. 2000) and luminosity-variability (Reichart et 
al. 2001). A particular relation between the lag and the variability had been strongly 
confirmed to prove both above luminosity indicators arc reliable (Schaefer ct al. 2001). 
On the other hand, other indicators based on GRB spectral features arc also constriictcd 
subsequently. They originate from cither the Ep — Eiao relation (Amati et al. 2002; Atteia 
2003), the Ep — Lp relation (e.g. Yonetoku et al. 2004) or the Ep — E^ relation (Ghirlanda 
et al. 2004b). The spectra and the light-curves are related to each other, via the spectral 
lag. 

Norris et al. (2004, 2005) found that wide pulse width is strongly correlated with 
spectral lag and these two parameters may be viewed as mutual surrogates in formulations 
for estimating GRB luminosity and total energy. Motivated by the above-mentioned 
developments, our first aim is to analyze the RSLs of long bursts in order to see what the 
distribution of them should be. Further purpose of this work is to search for the possible 
relations between some other parameters and them, and then interpret them in terms 
of physics. Data preparation is performed in § 2. In § 3, we are going to measure some 
typical physical variables. § 4 shows all the results. We are about to apply the RSLs to 
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observed data and try to reveal their physical explanations in § 5. We shall end with 
discussion and conclusion in § 6. 



2 DATA ANALYSIS OF PULSES 
2.1 Sample selection 

We use 64 ms count data selected from the current BATSE catalog for long bursts, called 
sample 1 including 36 sources. Note that we here only take into account these bursts with 
single pulse in the course of selection. The highly variable temporal structure observed 
in most bursts is deemed to be produced by internal shocked outflow, provided that the 
source emitting the relativistic flow is variable enough (e.g. Dermer & Mitman 1999; 
Katz 1994; Rees & Meszaros 1994; Piran, Shemi & Narayan 1993). In this case, the 
temporal structure generally reflects the activity of the "inner engine" that drives the 
bursts (Sari & Piran 1997). As a result of overlap, it is generally diflicult to determine 
how many pulses a complex bursts should comprise or to model the shape of these pulses 
(Norris et al. 1996; Lee at al. 2000). Fortunately, the observed peaks have almost one- 
to-one correlation with the activity of the emitting source, that is to say, each pulse 
is permissively assumed to be associated with a separate emission episode of one burst 
(Kobayashi et al. 1997; Kocevski et al. 2003). On the other hand, spectra parameters for 
distinct pulses within a burst are different from each other, which allows us to believe the 
spectral lags between these pulses will behave large difference (Hakkila & Giblin 2004; 
Ryde et al. 2005). Therefore we let our sample be composed of the relatively simple and 
bright bursts dominated by a single pulse event here rather than those dim or multi-peak 
ones for which we could accurately calculate the spectral lags. 

The method of selection is not automated program (e.g. Scargle 1998; Norris et al. 
2001; Quilligan et al. 2002) but directly experienced eyes, which in a certain degree could 
reduce any biases either from denoising techniques or from pulse identification algorithm 
itself (Ryde et al. 2003). Lee et al. (2000) has found the numbers of the pulses within a 
burst are usually different between energy bands. In principal, a bright-independent anal- 
ysis is required as the burst duration measurement needs (Bonnell et al. 1997), whereas 
the level of S/N should be reasonable and reliable. Based on these considerations, the 
criterions for our sample selection are now taken as follows: Tgo duration > 2.6 s; BATSE 
peak flux (50-300 kev) > 1.5 photons cm^^s^^; and peak coimt rate (> 25 kev) > 14000 
counts s~^. Next, we are going to process the data via background subtraction and 
denoising. 
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2.2 Background subtraction and denoising 

In general, the first step in data preparation is to select the appropriate background 
for substraction. To handle these data as a whole, an alternative mode of processing 
involving background subtraction along with denoising is presented here. For each source, 
we take the signal data cover its full range of the pulse as possible in order to ensure the 
contributions of all signals to lags are considered sufficiently. Prom the point of view of 
experience, data beyond this range are regarded as the fit of the background. However, 
for convenience, we here prefer disposing of the whole data involving pre- and post-pulse 
to processing separately. 

Considering the duration and background level of long bursts, we first smooth them 
with the DBS wavelet with the MATLAB software and then fit them with a pulse function 
plus a quadratic form, namely 

where the first expression on the right is a quite flexible function (see eq. (22), Kocevski 
et al. 2003) applied to describe pulse shapes and the hinder quadratic term represents a 
background that spans the whole data. The parameter tm is the time of the maximum 
flux, Fm , of the pulse and the quantities of r and d are next two indexes describing the rise 
and decay of pulse profiles respectively. Once the part of background is subtracted from 
the fitted data, the remainders are pure sinal data that isn't contaminated by background 
and noise in a certain error level. These signal data are just what we need to use for the 
analysis of spectra and light-curves. 

We divide our analyses into two portions in order to achieve different goals of calcu- 
lations for sample 1. One portion is to pretreat the light-curve data in energy channels 
1 and 3 (i.e., 25-55 kev and 110-320 kev). Here, we define these signal data belonging 
to the two channels as sample 2. Another is to combine the data from all four energy 
channels in order to study the characteristics of the "bolometric" light-curve profile, e.g. 
FWHM and asymmetry. To avoid confusion of the definitions, the summed signal data 
are called sample 3. Figure 1 displays two examples of the studied pulses (# 2387 and # 
6625). Panels (1) and (3) show the profiles in energy channels 1, 2 and 3 and their cor- 
responding fit curves identified by solid, dashed and dotted lines respectively. Panels (2) 
and (4) depict the the "bolometric" light-curve profiles and their fit curves symbolized 
by the smooth lines. 

3 PHYSICAL QUANTITIES MEASURE 

In this section wc focus our attention on measuring some quantities (RSLs, FWHM and 
asymmetry) of the single-peaked events from above signal data for 36 long GRB pulses. 
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Fig. 1 Light-curves of two of studied single-peaked events, GRB 930612 
2387, panels 1 and 2) with good S/N and GRB 980302 (# 6625, panels 3 and 4) 
with poor S/N. The best fit lines to observed light-curves in energy channels 1 
(25-55 Kev), 2 (55-110 Kev) and 3 (110-320 Kev) are denoted by smooth solid, 
dashed and dotted lines respectively (panels 1 and 3), and the "bolometric" 
light-curves combined from all four energy channels are shown with the best fit 
model (smooth curves) in panels 2 and 4 respectively. 

3.1 Relative spectral lag 

Applying sample 2 we cross-correlate energy bands between energy channels j and k with 
the following cross-correlation function (CCF) (Band 1997) 

CCF(r;.„..)= <"^-(^^-^(^ + -)> (, ^ fc) (2) 

where a^. =< vf >^/^, i [=1,2,3 or 4] represents different energy channel; r is the so- 
called spectral "lag" between any two of these channels; Vj and Vk stand for two time 
series in which they are respective light curves in two different energy bands. If the 
considered channels are appointed to he j = 3 and k — 1, the spectral lag can be thus 
written as T31, differing from those previous definitions of spectral lag (e.g. Norris et al. 
2000; Gupta et al. 2002), we otherwise define a quantity called RSL, namely 



= T3i/FWHM(i) 



(3) 
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where T31 represents the lag between energy bands 3 and 1; and FWHM(^i^ denotes the 
full width at half maximum of time profile in energy channel 1. The T31 is determined by 
the location of t where CCF peaks because the CCF curve on this occasion is smoothing 
and resembling gaussian shape near its peak. If the data points close to peak are not 
dense enough, we shall interpolate them within the range from one-side to another of the 
peak. One could find from this definition that Trei,3i is indeed a dimensionless quantity. 



3.2 FWHM &z Asymmetry 

For light-curve of a source, asymmetry and FWHM as the fundamental shape param- 
eters are especially required to be determined firstly since they can reflect properties 
of bursts themselves. The key factor of this issue may be that FWHM is associated 
with the energy of photons detected by observer (or the lorentz factors of ejecta with 
FWHM oc (see Qin et al. 2004)) and the asymmetry is largely influenced by co- 
moving pulse width if only the burst duration is not large enough, say, T90 < 1000s (Zhang 
& Qin 2005). Therefore investigations on asymmetry and FWHM should provide us a 
probe to detect their intrinsic parameters in bursts. 

We now commence measuring these parameters in a certain significance level (see § 
3.3). Note that unlike FWHM^i) (in eq [3]) FWHM here represents the full width at 
half maximum of the 'bolometric' light-curve profile. Asymmetry is defined by the ratio 
of the rise fraction (tr) of FWHM of pulse to the decay fraction (td). With a simple 
algorithm, we quickly measure the above quantities. Sometimes the sparse data may 
need to be interpolated with some points to improve the precision in calculation. For the 
purpose of analysis, samples 2 and 3 have been utilized distinguishingly. 



3.3 Error analysis 

Seen from eqs. (l)-(3), both fitted and derived parameters should be disturbed by errors 
propagated between these parameters directly or indirectly. It's necessary to illustrate 
the impact of background portion on observed variables has been eliminated from eq. (1) 
by background subtraction. Assuming the error of observed time (t) is zero at any one 
of data points, we find that the error of fiux should be caused by parameters F^, tm, r 
and d, i.e. 

In reverse, t can be expressed bs, t = t{x\,X2,Xi,Xi,X5) = f~^{F,Fm,tm,r,d), so the 
error of t derived by the fit to the observed data should be 

i=l ' i=l * 
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Suppose two ends of FWHM are ti and t2 respectively, then the errors of tr, td and 
FWHM can be written as 



a{tr) = [a^{h) + a^itm)]'/^ (6) 

<T{td) = [c7^{tm)+a\t2)]'/^ (7) 

a{FWHM) = [a^h) + a^{t2)]^/^ (8) 

where tr = tm — ti,td = i2 — im and FWHM = t2—ti. With the definition of asymmetry 
(= tr/td), its error can be gained by 

CF {asymmetry) = asymmetry x + ^ '^i'^d) ^2ji/2 ^g-j 

tr td 

The term v in eq. (2) actually represents the variable F in eq. (1) after background 
subtraction, which ensures that the error of CCF is caused by the fiuxes fitted to energy 
channels 1 and 3, that is to say 

a{CCF) = ii^fa^v,) + (^fa^vs)}'/' (10) 

where a'^{vi) and a'^{vs) are respectively determined by eq. (4) in energy channels 1 
and 3. In addition, the eq. (2) can be in turn replaced by rsi = T3-i{xi,X2,X3) = 
f-\CCF,vi,V3). Thus 

i—l * i—1 ^ 

It shows that the error of T31 comes from not only these fitted parameters but also the 
a{CCF) itself. The errors of Trei,3i can be deduced from eqs. (2) and (3) by the following 
propagation relation as eq. (9) 

..(.,.,30=.„u.x[(^)- + ( '™-'' )'lV- (12) 

Tai t W rlM(i^ 

where parameters FWHlv^i) and a{FWHM(^i)) in energy channel 1 can be measured 
with the same way as what is used in eq. (8) as well as eq. (5). 

In fact, for a given sample the rise time (f,.) and the decay time [td] have a linear 
correlation of tr oc C,td where the value of coefficient C, is estimated to be 0.3 — 0.5 for most 
bursts (Norris et al. 1996). The exact value of C, is uncertain because of the uncertainty 
in the parameters of the bursts in the sample. Considering the same reason, the relation 
between T31 and FWHM^^i), rgi « QM^FWHU^^^ '^'^ (Norris et al. 2005), is also required 
to be calibrated with a larger sample. Otherwise, we can't accurately determine the errors 
propagated from these connected parameters. So we assume they are independent. Under 
the circumstances, the errors caused by this un-correlation would become lager than that 
in the case of correlation. It's necessary to point out the assumption doesn't influence 
the credibility in our error analysis. 
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Fig. 2 Histogram of the distribution of RSLs for long bursts with a sample of 36 
bright pulses. The smooth curve is gaussian function fitted to the distribution, 
where the mean value is fj, — 0.102 and the standard deviation is cr = 0.045. 

4 RESULTS 

First of all, we display the RSLs distribution of long burst pulses (§4.1). To get the 
physical implications about it much detailed investigations have been given to one sub- 
sample (sample 4), which includes nine GRBs of long-lag and wide-pulse with simpler 
physics owing to more accurate measures (§4.2). 

4.1 The RSLs distribution 

Taking j = 3 and k = I and combining eqs. (2) and (3), we derive the quantities T31 and 
FWHM(i-f with samples 2, then Trei.ai- Plot of the RSLs distribution is illustrated in 
figure 2, from which we can find all the pulses hold positive Trei within the range from 
to 0.35 and they concentrate on the approximate value of 0.1. 

Moreover, we fit the distribution with gaussian function and get 'x^ /dof — 1.1 with 
= 0.97, which indicates that the distribution of RSLs is consistent with normal dis- 
tribution. However, the distributions of FWHM and spectral lags (or time intervals) are 
found to be lognormal instead of normal (see e.g. Mcbreen et al. 2003). In the following, 
we pay our particular attention to the properties of Trei,3i that would be studied in very 
detail as possible. 

4.2 The dependence of Trei,'ii on light-curve parameters 

To investigate whether Trei,3i is associated with some parameters of light-curve, we thus 
list the related parameters in table 1 and make plots of these relations as shown in figure 
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Table 1 Parameters gained by fit to sample 4 with the current model of this 
paper 



Trigger 


* 


T31 (sec) 


FWHM (sec) 


Asymmetry 


HR31 




tm (sec) 


1406 


0.126±0.017 


1.47±0.15 


9.64±0.30 


0.43±0.03 


0.85 


464.09±2.19 


4.47±0.05 


2387 


0.127±0.007 


2.43±0.04 


14.57±0.30 


1.38±0.07 


0.95 


738.89±2.03 


7.42±0.03 


2665 


0.166±0.014 


1.47±0.10 


4.56±3.40 


0.44±0.51 


0.45 


327.94±2.49 


3.22±0.26 


3257 


0.059±0.003 


1.61±0.07 


12.98±0.74 


0.27±0.02 


1.48 


462.97±2.02 


4.16±0.06 


6504 


0.105±0.006 


2.02±0.06 


9.78±1.13 


0.36±0.07 


1.49 


457.78±3.15 


4.39±0.11 


6625 


0.101±0.04 


1.55±0.19 


15.51±1.24 


0.57±0.08 


0.39 


338.18±2.61 


9.24±0.29 


7293 


0.092±0.009 


1.84±0.06 


12.11±2.66 


0.31±0.09 


1.74 


626.78±2.76 


5.32±0.15 


7588 


0.149±0.027 


1.14±0.22 


6.56±0.32 


0.56±0.05 


0.60 


449.81±2.63 


5.04±0.06 


7648 


0.192±0.036 


2.19±0.30 


11.94±0.60 


0.55±0.05 


2.12 


305±4.02 


6.89±0.13 



Notes: symbol * denotes the relative spectral lags calculated with equations (2) and (3), as 
shown in table 3. 



Table 2 Linear correlation analysis of figures 3 and 4. Coefficient and proba- 
bility are abbreviated to coef. and prob. respectively. 



Panels 






a 


b 


c 


d 


e 


f 


For figure 


3 


coef. 


0.08 


-0.48 


0.42 


-0.04 


-0.39 


0.03 






prob. 


0.84 


0.19 


0.25 


0.92 


0.29 


0.95 


For figure 


4 


coef. 


-0.40 


0.32 


-0.28 


-0.89 


-0.87 


-0.51 






prob. 


0.32 


0.43 


0.49 


0.001 


0.002 


0.16 



3. The hardness-ratio {HR31) here is defined as the ratio of photon counts of channel 3 
to channel 1 for pure signal data. The parameters Fm and tm as well as their errors are 
taken from the fit with eq. (1). 

For the sake of testing these relations, we take a general linear correlation analysis 
IjPress et al. 19921 and derive the correlation coefficients together with probabilities listed 
in table 2. In figure 3, unless the data point (# 7648) is removed from these panels (a, d 
and f), Trei,3i is evidently uncorrelated with T31, HR31 and t„j, even though it appears 
a tendency of reverse correlation between them. We find from the panels (b, c and e) 
that Trei,3i shows a positive correlation with asymmetry, while it is reversely correlated 
with FWHM and F,n- It demonstrates that the Trei.31 is very sensitive to FWHM and 
asymmetry, which suggests that it should be another parameter describing the shape of 
GRB pulses. 
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Fig. 3 Correlations between Trei,3i and other parameters of light-curves such as 
Tai, FWHM, asymmetry, HRsi, and tm- These relations will become more 
tight provided the source # 7648 is excluded from our sub-sample. 



5 APPLICATIONS AND EXPLANATIONS 

Encouraged by the special attributes of the Trej.si, we shall apply this quantity to explore 
some other potential correlations, so that we can have the opportunity to interpret it in 
terms of physical significance. 

5.1 Detecting relations between Trei,3i and other physical parameters 

First of all, we try to choose the appropriate variables meeting the aim of studies. Since 
how to determine the energy spectra as well as distances (or energies) is usually considered 
to be the key issue of understanding the burst phenomenon, we therefore select the 
relevant parameters such as a, /3, Ep, z, Lp and Fp to correlate with Trei,3i respectively 
(see table 3). 

It needs to clarify that the rcdshifts of all sources except # 7648 arc borrowed from 
their estimation by an empirical Ep — Lp relation instead of direct observations for the 
absence of their information of the afterglow. The reason for this selection is Ep — Lp 
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11 



Trigger 


'''rei, 31 


a" 




Ep 


z" 


Lip 




1 


2 


3 


4 


5 


6 


7 


8 


1406 


0.126±0.017 


-1.20±0.15 


-2.89±0.21 


79.43±21.39 


1.91±0.4 


22.8±10 


2.13 


2387 


0.127±0.007 


-0.03±0.09 


-2.46±0.03 


106.79±8.00 


3.76±0.4 


176±50 


3.58 


2665 


0.166±0.014 








1.19±0.11 


11±2 


1.91 


3257 


0.059±0.003 


-0.24±0.09 


-2.79±0.12 


169.2±16.32 


11.97±1.6 


1750±600 


2.62 


6504 


0.105±0.006 


0.59±0.26 


-2.78±0.17 


120.2±20.34 


1.67±0.07 


40±4 


2.25 


6625 


0.101±0.04 


-1.18±0.12 


-4.54±1.76 


53.8±9.31 


2.18±0.4 


15.8±6 


1.68 


7293 


0.092±0.009 


-0.15±0.12 


-2.81±0.11 


151.1±17.41 


8.48±3 


733±500 


2.73 


7588 


0.149±0.027 


-1.73±0.38 


-2.80±0.08 


12.07±17.9 


1.44±0.08 


7.79±2.2 


2.06 


7648 


0.192±0.036 


-0.78±0.26 


-2.43±0.21 


146.93±58.0 


0.43 " 


0.54±0.1 


1.53 



Notes.- Redshift (col. 6) and peak luminosity (col. 7) estimated by the Ep — Lp relation have 
been borrowed from Yonetoku et al. (2004) due to lack of the information about these sources 
except for trigger 7648 whose z and Lp (lO^'^ergs s~^) is offered by Galama et al. (1999) and 
Guidorzi et al. (2005) respectively. Note the unit of peak flux (Fp) is photons cm~^s~^. 
References.- a. Norris et al. 2005; b. Yonetoku et al. 2004; c. Galama et al. 1999; d. Guidorzi 
et al. 2005. 

relation looks considerably tighter and more reliable than the relations offered by the 
previous works |Yonetoku, et al., 2004| ), whereas the spectral parameters (a, /3 and Ep) 
for # 2665 are unavailable at present. 

All the correlations of T^ei.ai with these parameters are shown in figure 4. In the same 
way, we apply the linear correlation analysis to these relations and list the corresponding 
results in table 2. Figure 4 illustrates all parameters but /3 are reversely correlated with 
Trei.31- The high energy spectral index, f3, behaves an otherwise trend of positive corre- 
lation with Trei,3i- It has been proposed that Fp is often used as an effective indictor of 
distance to the GRB sources (Lee et al., 2000). The tight relation between Trei,3i and Fp 
in figure 4(f) shows that Trei.31 is expected to be a distance indictor. 

It is surprisedly found from panels (d) and (e) that T^ei.ai is strongly correlated with 
z and Lp respectively. Accordingly, our attempt in next section is to confirm if Trei,3i is 
suitable for an indictor of distance or peak luminosity. 

5.2 Redshift and luminosity indicator 

As Atteia (2005) points out, whether redshift indictors are good or not is determined by 
the degree of correlation between redshift and them which are generally combinations of 
GRB parameters with a small intrinsic scatter. To testify the validity of Trei.31 as the 
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1.0 - 




0.24 



, 1 , 1 , 1 1 1 , 1 1E50 I 1 1 1 1 1 1 1 1 1 1 I— I 1 1 1 1 1 1 1 1 1 

0.04 0.08 0.12 0.16 0.20 0.24 0.04 0.08 0.12 0.16 0.20 0.24 0.04 0.08 0.12 0.16 0.20 0.24 

'^rel.ai '''rel.ai '^rel,31 

Fig. 4 Correlations between Trei,3i and other parameters such as a, /3, Ep, z, 
Lp and Fp. The relations also have the tendency of becoming more tight than 
before once the source # 7648 is excluded from our sub-sample. 



redshift indicator, we contrast the observed data with the theoretical model in figure 
5. Note that the data of panels (a) and (b) are merely a replica of figure 4(d) and (e) 
correspondingly. Prom figure 5(a) and (b), the best fits to a linear model are written as 

log^;» 1.56-9.66T^ei,3i (13) 
log Lp « 55.44 - 23.07r„;,3i (14) 

with spearman rank-order correlation coefficients of -0.88 (p ^ 1.5 x 10^"^) for the former 
and -0.83 (p ^ 5 x 10^"^) for the latter. This indicates the relatively accurate connection 
between the RSLs and the redshift (or luminosity) does exist. Provided the Trei,3i is 
measured, using eqs. (13) and (14) one could precisely estimate redshifts and peak lu- 
minosities of those sources without the information of observed spectral lines. From this 
viewpoint, the quantity Trei,3i can be regarded as an ideal indicator of redshift and/or 
luminosity. 

Meanwhile, the RSLs might be utilized to constrain the cosmological parameters (say, 
f^TO) ^\ and Hq) in case redshift and luminosity are determined by eqs. (13) and (14) 
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Fig. 5 Calibration curves for the relative spectral lags, Trei,3i - Plots of T^ei.ai vs. 
redshift and peak luminosity can be used to calibrate redshift and/or luminosity 
indicators. The plots here can be fitted to yield eqs. (13) and (14) (marked with 
the straight red lines in panels a and b). 

simultaneously. Certainly, the realization of this purpose requires us to eliminate selection 
effects as possible as we can in advance, not only on observations but also on calculations. 

6 CONCLUSIONS AND DISCUSSIONS 

According to above investigations, we could derive the following conclusions: First, the 
distributions of RSLs are normal and concentrate on the value of 0.1 or so. Second, the 
RSLs are weakly correlated with FWHM, asymmetry, a, (3, Ep, Fp and While they 
are uncorrelated with lags^i, HR31 and tm- Finally, we find the Trei.31 is a useful redshift 
and peak luminosity estimator. 
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In view of internal shock model, GRBs generally consist of many pulses produced 
by multiple relativistic shells (or winds) followed by internal shocks due to collision of 
them as a central engine pumps energy into medium (e.g. Fenimore et al. 1993; Rees 
& Mcszaros 1994). Even so, our action on sample selection doesn't conflict with the 
idea of overlapping. For simplicity, we only choose single peaked bursts to construct our 
sample to make this issue easier to study. Furthermore, the manner of selection can also 
avoid the contamination of adjacent pulses by overlap, for which coidd inevitably produce 
additional errors, owing to selection effects. Under the assumption that these single wide 
pulses could be produced by the same mechanism (see Piran 1999), our motivation of 
this selection is advisable. In fact, McMahon et al. (2004) has pointed out the favorite 
mechanism for producing gamma-ray emission for such single pulse events could be not 
internal shocks but external shocks (see also Meszaros & Rees 1993; Sari & Piran 1996; 
Dermer et al. 1999, 2004). 

It is usually expected that the gamma-rays come from a relativistically expanding 
fireball surface with lorentz factor F > 100 (e.g. Lithwick et al. 2001). The bulk Lorentz 
factor increases linearly with radius as long as the fireball is not baryon loaded and 
not complicated by non-spherical expansion (Eichler et al. 2000) until F > 1000 (or 
R > IO^^cto) (Woods et al., 1995). With Power density spectrum method, Spada et al. 
(2000) found the curvature time together with the dynamic time will dominate over the 
radiative cooling time at a distance R < 5 x Given this case, based on Doppler 

effects the FWHM has been found to follow FWHM cx F^^ (qj^ 2004). Assuming 
T DC F^", from figure 3(b) we can deduce that the upper limit of ui is about 2 which is 
compatible with r cx F~^ gained by Shen et al. (2005). 

In this work we have shown that relative spectral lag can be used as an estimator 
of redshift and peak luminosity of long GRBs. We note that our conclusion is based on 
an analysis using nine sources. More accurate and robust results for the analysis would 
require a larger sample which includes sources with redshifts. Until now the expected 
indictor of short bursts hasn't been constructed yet, owing to lack of sufficient sources 
with the redshift, although redshifts of few sources have been measured by Swift and 
HETE II and reported recently (see, e.g. Bloom et al. 2005; Hjorth et al. 2005; Berger 
et al. 2005). The recent observations suggest that short bursts reside at cosmological 
distances, however, previous investigations showed the spatial population of short GRBs 
seems to accord with lower redshift sources (e.g. che et al. 1997; Magliocchetti et al. 
2003; Tanvir et al. 2005). Therefore, the research on redshift indictor of short bursts is 
excluded from this work and it still remains at the prenatal stage. 
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